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I .  ABSTRACT 


We  report  on  progress  of  work  supported  by  the  Air  Force 
Office  of  Scientific  Research  Grant  No.  AFOSR-80-0083  during  the 

\  I 

period  4/1/81  -  3/ 3 1/8 2."^ Three  problems  were  considered.  These 
were  finite  element  approximations  to  the  inhomogeneous  Navier- 
Stokes  equations,  for  potential  flows,  and  for  acoustic  eigenvalue 
problems.  In  all  cases  both  theoretical  error  estimates  and 
computer  codes  implementing  the  best  algorithm  were  developed. 

We  also  report  on  other  activities  sponsored  by  the  grant,  i.e. 
student  research  and  conference  talks. 
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II.  PROGRESS  REPORT 

1.  INHOMOGENEOUS  NAVIER  STOKES  EQUATIONS 

Finite  element  methods  for  stationary  viscous  incompressible 
flows  were  considered.  Specifically,  we  consider  the  problem  of 
finding  a  velocity  field  u  and  a  presure  field  p  which  satisfy 

V  •  u  =  g  in  ft 

(1)  vAu  -  u*Vu  +  Vp  =  f  in  0 

u  •  n  =  q  on  T 

2  3 

where  ft  is  a  bounded  region  in  3R  or  TR  with  boundary  r, 
and  where  g,  f ,  and  q  are  given  functions.  Previous  work  [1] 
on  the  topic  consider  the  case  g  =  0  and  q  =  0  only.  The  case 
q  £  0  is  important  since  usually  such  flows  are  driven  at  the 
boundaries,  e.g.  by  inflows.  Our  work  on  this  topic  included  the 
following  accomplishments: 

a)  Under  natural  hypothesis  on  f,  g,  and  q,  the  existence 
and  uniqueness  of  weak  solution  of  (1)  was  proven.  Adding 
stability  hypotheses  on  the  finite  element  spaces,  the  existence 
and  uniqueness  of  approximate  finite  element  solutions  was  also 
displayed. 

b)  Optimal  error  estimates  for  the  finite  element  discretiza 
tion  were  derived. 

c)  The  convergence  behavior  of  iterative  methods,  i.e.  the 
Newton, chord  and  simple  iteration  methods  for  the  solution  of  the 
discrete  nonlinear  algebraic  equations  resulting  from  the  finite 
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element  discretization  were  analyzed.  In  particular,  it  was  shown 
that  the  Newton  method  was  quadratically  convergent  when  one  is 
close  enough  to  the  exact  solution  of  the  discrete  equations,  and 
it  was  also  shown  that  a  simple  iteration  scheme  is  globaly  con¬ 
vergent  whenever  the. solution  of  (1)  is  unique. 

d)  Computer  programs  were  written  implementing  some  particular 
choices  of  finite  element  spaces.  These  programs  verified  the 
theoretical  results  and  also  served  to  illustrate  the  implementa¬ 
tion  of  the  algorithms  developed. 

These  results  were  reported  in  Paper  4  of  Section  III  below. 

A  crucial  step  in  developing  good  algorithms  for  the  ap¬ 
proximation  of  the  solution  of  (1)  is  choosing  finite  element  spaces 
which  satisfy  the  stability  hypothesis  alluded  to  in  (a)  above. 

We  choose  spaces  Vh  and  s'1  in  which  to  seek  an  approximation 
u  and  p  ,  respectively,  to  (1)  .  Then  the  stability  hypothesis 
takes  the  form 


£ 


*hdiv  vhdfl 


sup 

,h€V*  1^1 


R„ - >  Yhll»hH0  V*h  €  Sh. 


where  11*^  denotes  the  H1-Sobolev  norm.  The  crucial  question 
is  whether  or  not 


Yh  >  Y0  >  0 


uniformly  in  h,  i.e.  Yq  is  independent  of  h.  For  some  obvious 
choices  of  V*1  and  Yq  “  0  and  these  are  discarded.  For 

some  other  choices,  such  as  bilinear  velocities  and  constant 
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pressures,  Yq  J4  0  once  the  "checkerboard"  pressure  mode  is 
removed.  However,  there  is  some  evidence  that  y*1  depends  on 
h,  i.e.  y*1  “*  0.  This  results  in  possibly  bad  pressure  ap¬ 
proximations  which  must  be  filtered  in  order  to  obtain  useful 
pressures.  On  the  other  hand,  the  velocity  approximations  are 
optimally  accurate  without  any  filtering.  Computer  programs 
implementing  a  variety  of  low  order  elements  have  been  developed. 
These  programs  produce  accurate  velocity  approximations.  This 
work  is  reported  on  in  Papers  1,  2,  and  6  of  Section  III  below. 

2.  MIXED  FINITE  ELEMENT  METHODS  FOR  POTENTIAL  FLOWS 

This  work  involves  the  approximation  of  problems  of  the  type 


V  •  u  =  0  in  Cl 


u  =  V<P  in  Cl 


u  *  n  =  f  on  IV 


t  =  g  on  T 


where  T„  (1  L  =  T,  the  boundary  of  Cl,  f  and  g  are  given  func- 

N  U 

tions,  and  u  and  <P  are,  for  instance,  an  unknown  velocity 
field  and  potential  field.  The  simple  problem  (2)  is  equivalent 


At  *  o  in  Cl 


<j>  =  g  on  IV 


In  *  »  r0' 


However,  we  are  mtersted  in  discretizing  (1)  directly,  since  in 
more  general  settings,  problems  of  the  type  (2)  may  not  always  be 
recast  into  a  form  similar  to  (3)  .  Previous  work  [2]  on  finite 
element  methods  for  the  approximation  to  the  solution  of  (2) 
resulted  in  error  estimates  for  the  error  u  -  u*1,  uh  being 
the  discrete  solution,  in  the  norm 

Ihdl*  =  IIxIIq  +  Udiv  vii0 

where  |[»|L  denotes  the  L  -norm.  No  error  estimates  were 
obtainable  for  [Ju-u^JIq,  which  is  physically  of  much  greater 
interest  than  estimates  for  IIh-B^II**  The  thrust  of  our  work  was  to 
examine  conditions  under  which  optimally  accurate  approximations 
in  the  norm  11*  Uq  are  obtainable  for  the  problem  (2)  .  The  high¬ 
lights  of  this  work  are  the  following: 

(a)  Optimally  accurate  approximations  are  obtainable  whenever 
the  subspace  V*1  and  S*1  in  which  we  seek  our  approximate  solution 
u*1  and  <j>^,  respectively,  satisfy  the  inclusion  property 
Sh  =  divtV*1)  and  the  decomposition  property:  every  uh  ‘  €y>*  may 


be  written  in  the  fprra 


h  h  ,  h 
V  =  w  +  z 


where 


V  .  z  =  0 


0hlhshllo  <  Hdiv  wnLi 


where 


sup 

<i6hJ(Q) 

0*0 


J*f?aa 

fl 

mrm 


U  1- 

Again,  bounding  p  by  0n  >  0q  uniformly  in  h  is  crucial  to 
obtaining  optimally  accuate  approximations. 

(b)  An  example  of  a  pair  V11  and  Sh  satisfying  the  above 
properties. was  given. 

(c)  Computation  using  the  spaces  of  (b)  were  performed, 
verifying  the  optimal  theoretical  results.  Computations  using 
finite  element  spaces  which  violate  the  conditions  of  (a)  were 
also  performed,  and  non-optimal  or  divergent  approximations 
resulted.  This  work  was  reproted  on  in  Paper  2  of  Section  III. 


3.  EIGENVALUE  APPROXIMATIONS 

Work  was  completed  on  the  study  of  finite  element  approxima¬ 
tions  for  acoustic  eigenvalue  problems.  Here  we  consider  problems 
of  the  type 

A<t>  +  uj0  *  0  in  Q 

0  =  0  on  r 

recast  into  the  weak  form 

Ju-ydO  -  ^0divvdfl  =0  Vv  €  V 


(4) 


-£*div  udfl  =  V$  €  s 


for  suitable  Hilbert  spaces  V  and  S.  Error  estimates  for  finite 

element  approximations  of  (4)  were  performed.  Under  the  condi¬ 
tions  discussed  in  subsection  2  above,  these  estimates  were  optimal 


Computer  codes  implementing  stable  algorithms  were  developed, 
everywhere  verifying  the  theoretical  results.  This  work  was 
reported  in  Paper  1  of  Section  III. 
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